The r-matrix structure of the Euler-Calogero-Moser model 

E. Billey * J. Avan * O. Babelon * 

December 1993 



Abstract 

We construct the r-matrix for the generahzation of the Calogero-Moser system introduced by Gibbons 
and Hermsen. By reduction procedures we obtain the r-matrix for the 0{N) Euler-Calogero-Moser model 
and for the standard An Calogero-Moser model. 



PAR LPTHE 93-55 

*L.P.T.H.E. Universite Paris VI (CNRS UA 280), Box 126, Tour 16, 1"='' etage, 4 place Jussieu, F-75252 PARIS CEDEX 05 



1 Introduction 

The Calogero-Moser model seems to occupy a particular place in the world of integrable systems. On one 
hand, its complete solution is known ||^, ^ even at the quantum level, but on the other hand, the 

techniques which have now become standard — in particular the r-matrix formalism — did not seem to 
apply to it. 

The beginning of an answer to this problem was given in Q where the classical r-matrix was calculated 
and found to be a dynamical object, which is a non standard feature. This result was later extended to the 
elliptic case in |||, |^ . 

While the general theory of constant classical r-matrices is now well developped, this is not so for the 
dynamical ones and we still are in the phase of producing examples. Dynamical r-matrices appear when we 
consider higher Poisson structures ^ ^ . They also appear in non linear cr-models |ll| , or in the theory 
of Sine Gordon solitons |]l2| . Finally one should mention an interesting example of dynamical r-matrix for 
the hyperbolic Gaudin model |]l3| . 

In this paper we will consider the generalization of the Calogero-Moser model introduced by Gibbons 
and Hermsen iQ. This class of models has recently received some attention in relation with the evolution 
of energy levels with respect to the parameter t for a perturbated system of Hamiltonian H — Hq + tV 

Let us introduce a set of dynamical variables {qi,Pi)i=i..-N and {fij)ij=i-..N together with the Poisson 
brackets 



and the Hamiltonian 



{Pt,qj} = Sij 

{fijifkl} — Sjk fil ~ 5ii fkj 
1 ^ 1 ^ f- f- 



(1) 
(2) 

(3) 



In order to have a non degenerate Poisson bracket it is assumed that the {fij)i^j=i---N are restricted to a 
symplectic submanifold of (|2|). 

Notice furthermore that the quantities {fii)i=i---N Poisson commute with the Hamiltonian; we can there- 
fore reduce the flows to the surfaces {fa = constant)i=i...jv- These models were shown to be integrable 
precisely on these surfaces and our purpose is to compute their r-matrices. 

The plan of the paper is as follows : after a short review on the r-matrix formalism in section 2, we find 
in section 3 the Poisson structure of the Lax operator which implies the integrability property of the model 
on the manifolds {fa = constant)i=i...Ar. In section 4 we use the previous result to obtain the r-matrix 
structure for the 0{N) Euler-Calogero-Moser model introduced by Wojciechowski This is achieved by 
using the mean procedure introduced in [^, |2^, ^ . Finally in section 5 we again use the results of section 3 
to present a new construction of the r-matrix of the standard Calogero-Moser model computed for the first 
time in Q and obtained here by Hamiltonian reduction from the model (H). 



2 Dynamical r-matrices 

In the Lax representation of a dynamical system, the equations of motion arc written in the form 

L=[L,M] (4) 

where L and M take values in a Lie algebra Q with basis {^^^j. 

Eq.(^ implies that the quantities tr{L") are conserved. Liouville integrability requires their involution 
under the Poisson brackets. It was shown in |^ that this property {tr(L™), tr(L")} = is equivalent to 
the existence of an r-matrix. 

If we set 

Li = L(g)l 
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L2 = L 

this condition reads 

{ii,L2} = [ri2,Li]-[r2i,L2] (5) 

where 

ri2 = ^ r^^^Xfj, ® and r2i = ^ r^,.^,^ «) X^. 

In general the r-matrix is a dynamical object and possesses no special symmetry. Let us recall here some 
of its properties. 

1. The Lax operator L together with the r-matrix completely describe the dynamical system: if we choose 
iir(L") as a Hamiltonian, the equations of motion take the Lax form with M — tr2{ri2L2^^). 

2. The left-hand side of eq.(|^) being a Poisson bracket, the Jacobi identity implies a constraint on r. It 
is easy to check that it takes the form 

[Li, [ri2,ri3] + [ri2,r23] + [r32,ri3] + {i2,fi3} ~ {-^3,?'i2}] /gx 
-|- cyclic permutations = . 

It is interesting to note that the display of indices in (^) is not the one which appears in the usual 
Yang-Baxter equation. We recover it when the r-matrix is antisymmetric. 

3. Since the existence of an r-matrix is equivalent to the involution property of the eigenvalues of L, the 
form of eq.(^ is not modified by a conjugation of the Lax matrix. Indeed if L' = g~^Lg then 

{L\,L'2} = [r[2,L[]-[r'2„L'2] 

with 

r'i2 = 9i^92^ ~ {gi,L2}g^^ ^ ^["12,^2]^ 9x92 (7) 

where U12 = {9i,g2}9i^92^ ■ 

4. Finally we recall that the r-matrix of a given Lax matrix is not unique. In particular if we introduce 
a symmetric matrix T12 — T21 Q ® the change 

ri2 > ri2 + [ri2,i2] (8) 

does not affect the equation (||). 

3 The Poisson structure for the sl{N) model 

The Lax matrix of the system is 

i(A)^g(p.-f)e,.-f (-l-4)/,e„. (9) 

where {eij)ki = SikSji. 

The Poisson brackets of the elements of the Lax matrix (||) can be recast in the form 

^ f. - f 

{L^iX),L2{^i)} = [ri2(A,M),ii(A)] - [r2i(M, A), L2(^)] - J2 /'"_ ^'2 ^ 

where 

C ^ 1 

ri2(A,/x) = - ^ «) eji (11) 

^~ ^- --z^ Qi ^ 1] 
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and C is the Casimir element of sl{N) 

N 



C — ^ ^ Cij (S) Gji- 



Some comments are in order: 



1. Eq.(|To|) holds for the non reduced dynamical system (|^) and it is known that trL"{X) are not in 
involution for this system. This fact is responsible for the additional term in eq.(lO) compared to 
eq. (|). Indeed we have 



^ f---f- 



2. The constraints fa generate on L{X) a conjugation by a diagonal matrix. Hence trL"'{X) Poisson 
commute with fa. It follows that one can compute their reduced Poisson bracket on the manifold 
ifii = constant)i=i...Ar. This amounts to setting fa = constant in eq. (|lO|) and therefore trL^{\) are in 
involution. 

3. Since the quantities trU^iX) commute on the reduced phase space, L[\) has an r-matrix. However its 
explicit computation requires some care since L[X) is not a function on the reduced phase space. 

In the following two sections we give two examples of reductions where we can explicitely obtain the r-matrix, 
starting from the initial structure (|l0|). 



4 The r-matrix for the 0{N) model 

We now consider the following model, introduced by Wojciechowski in The dynamical variables are 

i<liyPi)i=i---N and antisymmetric {hij = — /iji)ij=i- -JV- The Poisson bracket is now: 

= % (12) 
{hij,hki} = - (Sii hjk + Ski hij +6jk hu + Sij hki) ■ (13) 

The Hamiltonian of the system reads 



As shown in |Tq | the equations of motion can be written in the form 

L(A) = [L(A),M(A)] 

where the Lax pair (L, M) is given by 



^ ' 1 1 



^(A) = - h - I hij eij (14) 



N 



and the Hamiltonian is H = ^ J ■^jxtrL{\)'^ . 

In ||l8| this system was proved to have Poisson commuting Hamiltonians trU^{\). It follows that this 
model admits an r-matrix, which we now calculate. 
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First of all we remark that the Lax matrix of the 0{N) model is obtained from the Lax matrix of the 
sl{N) model as 

LOW(A) = i (I-(t)L^'W(A) (16) 
where a is the following involutive automorphism of the loop algebra: 

a : A"ey ^ -(-A)"ej,. (17) 

In this average procedure 



1 

2 

and the Poisson structure (§) of the sl{N) model becomes the structure (^ 
It follows immediately that 



= 7;if., - fjd (18) 



^ i [It Qj) ' 

We now remark that 

^ f - f 

{l-a®l-l®a + a®a)Y^ ^ ® eji 0. (20) 

Moreover since 

(T(8)cr ri2(A,^) = -ri2(A, ^) (21) 
we finally get an explicit r-matrix structure for L'^^^\X): 

{L?(-)(A),L,°(-)(;.)} = (A,^),L?(-)(A)] - [r^^ (m, A), L,°(-)(m)] (22) 

with 



1 
2 

or explicitely 



rir^(A,A^) - -(1 + a® l)r^r^(A,M) (23) 



0(N) 

N 



- iE(^+A^)^'.^^- (24) 

This is an example of application of the well known mean procedure ||l^, |2^, |2^ . 

It is interesting to see how the Jacobi identity (^) is fulfilled for this dynamical r-matrix (p^). We find 
that 

1 ^ 1 

['"I2,ri3] + [ri2,r23] + [r32,ri3] = ^ X! 7 — Z — ^2^»J ® K^*:*' + '^:''*) ® ~ '^Jj) ~ ("^^ ~ "^w) ® ('^^ + ■ 
It is then easy to check that 

[''12, 7'13] + [''12, ?'23] + [''32, ''13] + {-^2, ''is} ~ [L^, ''12} = 0. (25) 

The remarkable feature of this equation is the vanishing of its right-hand side. This is in contrast with other 
cases [^ ^ where eg. (p5|) has a non vanishing right-hand side, albeit of a special form. 



4 



5 A new construction of the r-matrix for the Calogero model 

We now show how the dynamical r-niatrix for the usual Calogero model ^ stems from eq. ( p^ through 
a procedure of Hamiltonian reduction. 

Let us restrict oneselves to the symplectic manifold 

fti = 6: Vj (26) 

with = 0, {ili,r]j} = 0, {^i,T]j} = —Sij. The Poisson brackets of fij are indeed given by (||). 

On this manifold we have a symplectic action of an Abelian Lie group 

i^ K^m- (27) 

The Hamiltonian ^ is invariant under this action and we can apply the method of Hamiltonian reduction. 
The infinitesimal generator of this action is 

N 

He = ^ Ci/ii- 

1=1 

The reduction is performed by first fixing the momentum. We choose: 

The isotropy group Ga of a is the whole group itself since it is Abelian and we still have to quotient out the 
action of this group. At the end of the procedure, the 27V degrees of freedom {£,i,ili)i=i---N are eliminated, 
leaving as reduced Hamiltonian (|^) the Hamiltonian of the usual Calogero-Moser model 



1 ^ 1 " 1 

''^^'^2&^-2"^g(^-^■ 

In order to perform the reduction at the level of the Lax matrix we remark that the group acts on L as 
conjugation by the matrix diag(Ai)i=i...jv- We now observe that the matrix 

icai = diag(^,^^) L diag(^j) 
is invariant under Ga and we can therefore compute the Poisson brackets of its matrix elements safely. 

^ ^ 1 1- 



i=l i,j = i 



The r-matrix of this L matrix is now obtained from eq.(^0|) and (11) by applying the conjugation formula (|^). 
When setting fa — a the extra term in eq.(lO) vanishes, leaving us with the r-matrix 

1 ^ ^ / 1 i\ 
r^f ( A, ^) = ri2 ( A, ^) eii » e^i -I- ] ® eji 

which gives 

r'^2^{\ H) = — r en (g) eii - ( — ^- h — ) Cy (g) ej^ + ( — ^ ) e^i ® e^j. (29) 

This r-matrix does not coincide with the r-matrix obtained in |§, 0] , however it is easily related to it by 
a transformation (^) with T12 = ^ Tlii^i ® ^a- 
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6 Conclusion 

The class of models introduced by Gibbons and Hermsen turns out to be a particularly interesting 
generalization of the Calogero model. In particular their r-matrix structure is rather elegant and is general 
enough to allow the construction of a host of new models by various reduction procedures. Generalization 
to trigonometric and elliptic potentials is possible and will be developped further. 
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